A NECESSARY SOLUTION 
CONDITION FOR SUDOKU 



THOMAS FISCHER 

Abstract. We develop a new discrete mathematical model which 
includes the classical Sudoku puzzle, Latin Squares and gerechte 
designs. This problem is described by integer equations and a spe- 
cial type of inequality constraint. We consider solutions of this 
generalized problem and derive a necessary condition on these so- 
lutions. The results are illustrated with examples. 



1. Introduction 

A Sudoku is a square consisting of a 9x9 grid which is partly pre- 
populated by numbers between 1 and 9 called the givens. The problem 
consists of finding numbers between 1 and 9 for all unpopulated cells, 
such that each row, each column and each block consists of exactly the 
numbers 1, . . . , 9. The blocks of a Sudoku partition the Sudoku square 
into subsquares of size 3x3. Each Sudoku consists of 9 rows, 9 columns 
and 9 blocks. 

Example 1.1. A lot of examples for Sudoku are spread over the math- 
ematical literature and you can also find some in your daily newspaper 
or spread over the internet. We refer here to a typical Sudoku square 
of Delahaye [3J depicted in Fig. [H 

In Section [2] we introduce a mathematical model based on an integer 
constraint system, which we call the generalized Sudoku problem. In 
Section [3] the Sudoku puzzle is reformulated as a special case of this 
problem. 

The mathematical model depends on an integer n > 2 and is more 
general than a Sudoku puzzle. Sudoku problems of arbitrary size are 
covered by this model. The present approach also covers Latin Squares 
(see Denes and Keedwell [I]) and "gerechte Anordnungen" , introduced 
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Figure 1 . A typical Sudoku 



by Behrens p]. The latter is usually termed as gerechte designs (see 
Vaughan [T5] ) in the literature. 

The main result in Section [5] states a necessary condition on the so- 
lutions of the generalized Sudoku problem. The proof of this result is 
based on algebraic properties and does not involve a separation theo- 
rem. 

The use of necessary conditions is a powerful tool in the investiga- 
tion of nonlinear optimization problems (see Fiacco and McCormick 
[5] or Luenberger [9]). Beside optimization theory the Gaussian elim- 
ination process for solving linear systems of equations starts with the 
assumption of an existing solution and derives necessary conditions on 
this solution. In the theory of discrete optimization usually other tech- 
niques like relaxation or cutting planes are used (see Nemhauser and 
Wolsey [10]). 

There exist several papers dealing with algorithms for solving Su- 
doku. Classic approaches are brute force methods or paper-and-pencil 
methods (Crook [2]). It is also possible to use branch-and-cut algo- 
rithms (see Kaibel and Koch [7]) or the Algorithm X implementing 
dancing links described by Knuth [8]. 

The present paper describes a Sudoku puzzle as an integer program- 
ming problem. Similar approaches can be found in the literature (see 
Fontana |6]). To the best of my knowledge there do not exist results 
or methods based on necessary solution conditions, so far. 
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Finally, we collect some basic terms and notations. Let Z denote the 
set of integers and let Z™ denote the n-times cartesian product. The 
elements of a set are called distinct if they are pairwise distinct, i.e., 
each two elements of the set are not equal. The vectors respectively 
1 denote the zero respectively one vector, consisting only of zeros re- 
spectively ones in each component. The number of components is often 
indicated by an index. Each vector is considered to be a column vector. 
U denotes the identity matrix. The transpose of a vector or a matrix 
is indicated by a superscript T. We allow explicitly the use of vectors 
respectively matrices of dimension zero, i.e., without components re- 
spectively entries. In this case we call the matrix the empty matrix. 
The sign function is denoted by sgn and is generalized in Definition 14. II 
to vectors. We consider the sum over an empty index set to be zero. 



2. The Mathematical Model 
Let n be an integer with n > 1. We define the sum 

n-1 
i=l 

and define a matrix A(n) with s(n) rows and n columns inductively. 
For n = 1, let A(l) denote the empty matrix, i.e., a matrix without 
entries. Assume the matrix A(n — 1) had been defined with s(n — 1) 
rows and n — 1 columns. Then we set 





\ 


3-n-l 




, Os(n-l) 


A(n-1) ) 



Example 2.1. The matrix A{ri) for n = 9 is depicted in Fig. 0. 

We state two elementary properties of the matrix A(n). 

1. The vector of ones l n = (1, . . . , 1) T £ Z n is mapped to the zero- 
vector by A(n), i.e., A(n)l n = S („). 

2. The matrix A(n) has rank n — 1. 

We extend the matrix A(n) to a matrix A with n ■ s(n) rows and n 2 
columns. The matrix A consists in the "main diagonal" of n matrices 
A(n) and the remaining values are set to zero. The matrix A depends 
on the value n, but we do not state this dependence explicitly. 
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Given the set {1, . . . , n 2 } C Z, let n be any permutation on this set, 
i.e., 

n:{l,...,n 2 }^{l,...,n 2 } 

be a permutation. We extend the notion of permutation to the matrix 
A, i.e., we define n(A) = (a 71- ^, . . . , a 71 " ( n )), where a- 7 denotes the 
column of A for j = 1, . . . , n 2 . Analoguously we define the term 
ix(x) = (2^-1(1), . . . ,x n -i(n2)) T , for any x = (x u . . .,x n 2) T G Z™ 2 . We 
define the matrix A n = ir(A), i.e., we interchange the columns of A 
according to the permutation ir. 

Definition 2.2. Let s > 1. For any point y = (y 1 , . . . ,y s ) T G Z s we 
write y <> if each component of y is nonzero, i.e., if yi ^ for 
i = l,...,s. 

This definition should not be confused with the expression y ^ 0, 
where only one component of y has to be nonzero. 

Now we are in position to state the generalized Sudoku problem. 
Given is n > 2, some permutations tti, n 2 , n 3 on {1, . . . , n 2 }, some < 
k < n 2 , an index set {ii, . . . , i^} C {1, . . . , n 2 } and givens g ix , . . . , g ik G 
Z with 1 < g it < n for I — 1, . . . , k. 

The generalized Sudoku problem consists in finding x — (xi, . . . , x n 2) T 
in Z n , such that 1 < x-i < n for i = 1, . . . , n 2 , 

A ni x <> 0, 

A^z <> 0, 

and Xi t — g% t for I = 1, . . . , fc. 

We will investigate the relation between this model and the Sudoku 
puzzle in detail in the next section. 
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Figure 2. The matrix A(9) 



Due to the prescribed givens some of the components of x are known 
in advance. Therefore it is possible to reduce the size of this problem 
from n to n — k variables. This is advantageous for numerical purposes, 
but may interfere theoretical properties. We decided not to eliminate 
these variables in our considerations and maintain a distinction between 
equality and inequality constraints. 

This general approach includes gerechte designs (see Vaughan [13]). 
In particular the case ir^ = tt2 is allowed and in this case the generalized 
Sudoku problem describes Latin Squares (see Denes and Keedwell [1]). 

The generalized Sudoku problem can be classified as a nonlinear inte- 
ger problem. This problem even possesses nonconvex constraints. This 
contrasts with famous integer optimization problems like the traveling 
salesman problem and the knapsack problem, which usually are mod- 
eled as linear integer problems (see Papadimitriou and Steiglitz [TT]). 
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3. Modeling Sudoku 



In this section we show, that a classical Sudoku puzzle can be re- 
fomulated as a generalized Sudoku problem. We state an elementary 
result for later use without proof. 

Lemma 3.1. Let n > 2 and Xi,...,x n be integers. The following 
statements are equivalent: 

(i) xi, . . . , x n are distinct and 1 < Xi < n for % — 1, . . . , n. 
(li) {x 1 , ...,x n } = {l,...,n}. 

Lemma 3.2. Let n > 2 and x — (x 1 , . . . , x n ) T G Z n . Then A{n)x <> 
if and only if the components X\ , . . . , x n of x are distinct. 

Proof. We show the claim by induction on n. The claim is true for 
n = 2, since A(2)x = x\ — X2 and this expression is nonzero if and only 
if x\ and x% are distinct. 

Assume the claim holds for n — 1 and consider x = (xi, . . . ,x n ) T . 
By induction assumption, A{n — l)(x 2 , . . . , x n ) T <> if and only if 
distinct. Obviously distinct if and only if 

%i — x i 7^ for % = 2, . . . , n and distinct. The condition 

%i — Xi 7^ for i = 2, . . . , n can be expressed as 



We combine the Lemma 13.11 and 13.21 and obtain: 

Lemma 3.3. Let n > 2 and let x = (x\, . . . , x n ) T G Z n . The following 
statements are equivalent: 

(i) A(n)x <> and 1 < X{ < n for i = 1, . . . , n. 

(ii) {x 1 , . . .,x n } = {1, . . . ,n}. 

We introduce the notion of an x-tableau, where the components 

2 

xi, . . . ,x n 2 of points x G Z n are arranged row- wise from the left in 
an n x n square. 



(1 



n— 1 




<> 0. 



By definition of A(n) the claim is proved. 



□ 
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In the sequel we do not distinguish between the x-tableau consisting of 
the components of x and an x-tableau consisting of the corresponding 
indices. 

Let 7Ti be the identical permutation on {1, . . . , n 2 }, i.e., 7Ti(i) = % for 
% = 1, . . . , n 2 . If we arrange a point x = (xi, . . . , x n 2) T £ Z n , such that 
1 < Xi < n for i = l,..., n 2 , in an x-tableau, then A ni x <> holds if 
and only if the values in each row of the x-tableau consist of 1, ... ,n. 
This follows from Lemma 13.31 and the definition of A ni . 

From the next lemma we need part (i) in this section and part (ii) 
in a later section. 



Lemma 3.4. For any permutation ir on {1, . . . , n 2 } holds: 

(i) A n x = Atc~ 1 (x) for each x £ Z™ 2 and 

(ii) Al\ = ir(A T \) for each X £ {-1, +l} n < n \ 

Proof, (i) Let x = (x\, . . . , x n 2) T £ Z™ 2 , then 



A n x = ^^a n = rfXirtf) = An 1 ( 

3=1 3=1 



X). 



(ii) Let A £ {-l,+l} n < n \ then 





(a^X\ 






AlX = 




= 7T 










\a n2 Xj 



tt(A t X). 



□ 



s 
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Let 7T2 be the permutation, which maps the rows of the x-tableau to 
the columns of the x-tableau, i.e., 

7T 2 (1) = 1 

7r 2 (2) = n + l 
tt 2 (3) = 2 -n + l 

7r 2 (n) = {n — 1) • n + 1 
vr 2 (n+l) = 2 

Using Lemma I3.4[ if we arrange a point x = (xi, . . . , x n i) T G Z" 2 , 
such that 1 < Xi < n for z = 1, . . . , rz 2 , in an x-tableau, then A W2 x = 
Air^ 1 ^) <> holds if and only if the values in each column of the 
x-tableau consist of 1, . . . , n. 

If n is a square number the x-tableau can be covered by n non- 
overlapping subsquares of size \fn x \fn and we number these sub- 
squares row-wise from the left. Let 7r 3 be the permutation, which maps 
the z-th row of the x-tableau to the z'-th subsquare for i = 1, 



TTsjv^+l) = n+l 

n 3 (2-y/n) = n + ^/n 

7r 3 (n) = (y/n - 1) ■ n + y/n 



Using Lemma 13.41 if we arrange a point x = (x±, . . . , x n 2f G Z n , 
such that 1 < Xi < n for i — 1, . . . ,n 2 , in an x-tableau, then A n3 x = 
A-k^ 1 (x) <> holds if and only if the values in each y/nx y/n subsquare 
of the x-tableau consist of 1, . . . , n. 

If n = 9 and ttx, 7r 2 , are defined as before the generalized Sudoku 
problem is equivalent to the classical Sudoku problem. Each solution 
of the original Sudoku problem solves the generalized Sudoku problem 
and vice-versa. 

Also the Sudoku puzzle can be modeled as a linear integer problem 
(see Kaibel and Koch [7] and Provan [12]). Kaibel and Koch [7] pro- 
posed a model where the Sudoku problem is described by 0-1- variables 
and gave a short overview on important research results without stating 
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results explicitly. They also described how to apply existing software 
packages to their problem successfully. The same model had been con- 
sidered by Provan [12], who proved a unicity result. 

Another model is to formulate a Sudoku puzzle as an exact cover 
problem. The matrix describing this problem is the transpose of the 
matrix considered by Kaibel and Koch [7] and Provan [12] . 

4. Preliminaries 

In this section we start with a generalized sign function based on the 
classical sign function, also denoted by sgn. 

Definition 4.1. Let s > 1. For any point y = (yi, ... , y s ) T G If with 
y <> we define the generalized sign function sgn : If — > 1/ by 



'sgn(yi)\ 
K sgn{y 



■ 



sgn(y) = 

We return to the matrix A{n). 

Lemma 4.2. Let n > 2 and ) T G Z n with A(n)x <> 0. 

The point A T {n)sgn(A{n)x) G TP 1 consists of the components 

i— 1 n 

— sgn(xj — x,j) + sgn{xi — Xj) 

3=1 j=i+l 

for i — 1, . . . , n. 

Proof. We prove this lemma by induction on n. For n = 2 and x = 
(xi,x 2 ) G Z 2 with xi 7^ x 2 , 



A T (2)sgn(A(2)x) 



-l) S9H{Xl - X2) = (-sgn(xi -t)) 



On the other side 

2 

— sgn(xj — xi) + sgn{x\ — Xj) = sgn(x\ — x 2 ) 

3=1 3=2 

and 

1 2 

— sgn(xj — x 2 ) + sgn{x 2 — Xj) = —sgn{x\ — x 2 ), 
j=i j=3 

which proves the claim for n = 2. 
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Assume the claim is true for n — 1. Let x — (x±, . . . , x n ) T G T/ 1 with 
<> 0. The induction assumption reads as 



'x 2 



A T (n- l)sgn{A{n- 1) : ) 



( 



- sgn(xj - x 3 ) 

j=2 



n v 

J'=3 
n 

X) sgn(x 3 - Xj) 

j=4 



n-2 



sgn(xj - Xn-i) + Yl sgn(x n -! - xj) 



\ 



i=2 

n-1 

- E sgn(xj - x n ) 

3=2 



j=n 



By definition of A(n) 



( sgn(xi - x 2 ) \ 



sgn(A{n)x) 



\ 



sgn(x l - x n ) 

(X2 
x n / ) 

with n — 1 components sgn(x\ — Xj) for i = 2, . . . ,n and a vector 
sgn(A(n — l)(x2, ■ ■ ■ , x n ) T ) with s(n — 1) components. Again by defi- 
nition of A(n) 



A T (n) 





0f(n-l) 




A T (n — 1) 



We obtain 



A T (n)s^n(A(n)x) 



X) spn(a:i - Xj) 



sgn(xi — rc 2 



+ A T (n - l)s^n(A(n - 1) 



^2 
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E sgn(x 1 - Xj) 

3=1 



-sgn(x 1 - x 2 ) 
-sgn(x 1 - x 3 ) 

-sgn(xi - x n -i) 
-sgn(xi - x n ) 
( 



+ E sgn(x 2 Xj) 

3=3 



E sgn(xj -x 3 ) + E sgn(x 3 - x s ) 



J'=2 



n-2 



J=4 



E s^n(xj - x n -i) + E sgn(x n -! - Xj) 



3=2 
n-1 

E sgn(xj - x n ) 
i=2 



j=n 



n 

E S5m(a;i - X,) \ 



J'=2 



—sgn(xi — x 2 ) + E sgn(x 2 - Xj) 

2 ra 

- E sgn(xj -x 3 ) + E sgn(x 3 - Xj) 



n-2 



j=4 



E s</n(xj - x n _i) + E s#n(x n _i - a^) 



n-1 

" E - X n) 

3=1 

and this shows the claim. 



j=n 



\ 



□ 



Lemma 4.3. Let x = (xi, . . . , x n ) T G Z n wift {xi, . . . , x n } = {1, . . . , n}. 
JTien 

i— 1 n 

— sgn(xj — Xi) + sgn(x,i — Xj) = 2xi — (n + 1) 

jr' = l j=i+l 

/or i = 1, . . .,n. 

Proof. We prove this lemma by induction on n. For n = 1 both sides 
of the equation equal zero. This proves the claim for n = 1. 

Assume n> 2 and the claim is true for n— 1, i.e., for each (x 1 , . . . , x n _i) T G 
Z n_1 with {xi, . . . , x n _i} = {1, . . . , n — 1} we have 

i—l n— 1 

— sgn(xj — x^ + sgn(xi — Xj) = 2xi — n 

j=l j=i+l 

for i — 1, . . . , n — 1. 
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Let x G Z n with {x±, . . . , x n } = {1, . . . , n). Select the index 1 < 
i < n, such that x io = n and define x' G Z n_1 by 

X (^1) ■ • • j %io— 1) ^io+lj • • • j 2-n) • 

Let i G {l,...,n} and we distinguish 3 cases. 
Case 1: i < i - Then 

i— 1 n 

- ^ sprafo- - Xi) + ^2 Sf jn(xi - Xj) 

j=l j=i+l 
i-1 

3=1 

iq — 1 n 

+ ~ x i) + - ^o) + 5^ s 9 n ( x i ~ x j) 

j=i+l j =io+l 

i— 1 »i— 1 

= - XI s # n (4 - x i) + X s^n(x- - a£) - 1 

j=l j=i+l 

=2x'i -n-1 
=2xi - (n + 1). 

Case 2: i > i . Then we set / = z — 1 and obtain 

j— 1 n 

— sgn(xj — Xi) + sgn{xi — Xj) 

3=1 j=i+l 

io— 1 i— 1 

= - ^2 s 9 n ( x j ~ x i) - sgn(x io - x^ - ^2 sgn(xj - xi) 

3=1 j=io+l 



+ X sgn(xi - xj) 



3=i+l 

i-1 n-1 

= ~Y1 s # n K- ~ x 'i-i) + Yl s 9 n ( x 'i~i - x 'j) - 1 

3=1 j=i 
l-l n-1 

= ~Y1 s 9 n ( x 'j ~ x i)+ Yl sgn ( x 'i ~ ~ 1 

3=1 3=1+1 

=2x\ 

=2^-1 - « " 1 
=2xi -(n + 1). 
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Case 3: i = %q. Then Xi = Xi = n and 



i-l 



- 2J sgn(xj - Xi) + 22 s 9 n ( x i - 



X 3, 



i—1 n 

= — sgn(xj — n) + sgn{n — Xj) 
3=1 j=i+i 

=(i - 1) + (n - i) 

=n — 1 

=2n- (n+1) 
=2x, — (n + 1), 

and this completes the induction proof. □ 

We combine the Lemma 14.21 and 14.31 with Lemma 13.11 and 13.21 and 
obtain 

Lemma 4.4. Let n > 2 and ) GZ" components 

1 < Xi < n for i = 1, . . . , n and /e£ <> 0. Then 

x = ^(A 7 \n)sgn{A{n)x) + {n+ l)l n ). 

Proof. Because of Lemma 13.11 and 13.21 we are able to combine the 
Lemma 14.21 and 14.31 and we obtain 

A T (n)sgn(A(n)x) = 2x - (n + l)l n 

and the rest is elementary arithmetic. □ 

Example 4.5. We consider the case n = 9. The matrix A(9) is de- 
picted in Fig. [H and we choose the point x = (2, 8, 1, 5, 9, 4, 6, 3, 7) T . 
Then sgn(A(9)x) 

= (-1, 1, -1, -1, -1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 

-1, -1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -If 
and finally 

^(A T (9)sgn(A(9)x) + 10-l 9 ) 
= -(-6, 6, -8, 0, 8, -2, 2, -4, 4) T + 5 • 1 9 = (2, 8, 1, 5, 9, 4, 6, 3, 7) T = x. 
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5. The Necessary Condition 



In this section we extend the condition on the matrix A(n) in Section 
H] to the "large" matrices A n (see remarks prior to Definition 12. 2j) for 
any permutation ir on {1, . . . , n 2 }. 

Theorem 5.1. Let n > 2, let ir be any permutation on {l,...,n 2 } 
and let x = (xi, . . . , x n -i) T G Z™ 2 with components 1 < X{ < n for 
i — 1, . . . , n 2 and A n x <> 0. Then 



x 



-(^sgn{A K x) + (n + 1)1^). 



Proof. In the first step we prove the theorem for A n = A, i.e., n = id, 
the identity mapping. Let ) GZ™ with components 

1 < Xi < n for i - 1 ,J 



l,...,n and Ax <> 0. Using the diagonal 



structure of A and Lemma I4.4[ we obtain 

sgn(A(n 



and 



sgn(Ax) 



X(n— l)-n+l 



sgn(A(n) 



\ 



( 



A T sgn(Ax) 



A 7 \n)sgn(A(n 




X n'2 



x (n-l)-n+l 



A T (n)sgn(A(n) 



V 



Xri.2 




\ 



-(n+l)l n 




2x - (n + l)l n 2. 
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This shows the statement for A n = A. 

Now let 7r be an arbitrary permutation on {l,...,n 2 }. Let x = 
{x\, . . . , x n 2) T G Z n with components 1 < Xi < n for i = 1, . . . ,n 2 and 
A w x <> 0. 

Using Lemma [3.41 (i) An~ 1 (x) = A n x <> and we apply the result 
of the first part to 7r~ 1 (x). Now using Lemma [3.41 {it) we obtain 

A^sgn{A^x) = ir(A T sgnA-K^ 1 (x)) 

= n(2n- l (x) - (ra + l)l n a) 
= 2x — (n + l)l n 2 
and this completes the proof. □ 
This result can be applied to the generalized Sudoku problem. 

2 

Theorem 5.2. Let n > 2 and let x G Z n be a solution of the gener- 
alized Sudoku problem. Then 

x = 2^r S 9 n ( A ^r x ) + + 1)1„2). 
for each r = 1,2, 3. 

Proof. This follows immediately from Theorem 15.11 and the definition 
of the generalized Sudoku problem. □ 

Example 5.3. We consider the case n = 3, which implies s(n) = 3 
and consider a solution i6Z 3 which is depicted in this square: 



The matrix A. 









2 


1 


3 






3 


2 


1 


1 


3 


2 


in 


this 


case 


looks like 






rl 


-1 

















1 


1 








-1 

















1 





-1 


























1 




1 

















1 





-1 




















1 


-1 











1 


ii 


- 


1 














1 














-1 


.0 











1 








-1 



sgn^A^x) = (1, —1, —1, 1,-1 — 1,1,1, — 1) T and finally 

\{Al z sgn{A^x)+A - 1 32 ) 
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-(0,-2,2,2,0,-2,-2,2 



0) T + 2-l 9 = (2,l,3,3,2,l,l,3,2) T 




Each solution of the generalized Sudoku problem coincides with the 
givens. This can be used in a reformulation of Theorem 15721 The brack- 
ets with an index in the subsequent corollary denote the corresponding 
component of a vector. 

2 

Corollary 5.4. Let n > 2 and let x e Z n be a solution of the gener- 
alized Sudoku problem. Then 
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for each I = 1, . . . , k and each r = 1, 2, 3. 
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